CLAIMS 

What is claimed is: 



1 . A method for color matching a first image and a second image, wherein a first region of 
the first image and a second region of the second image overlap, the method comprising: 

generating a first histogram of the first region; 

generating a second histogram of the second region; 

determining corresponding pixel values firom the first and the second histograms; 

determining at least one parameter of an optoelectronic conversion fiiiiction 
(OECF) that best matches the corresponding pixel values; and 

color matching the second image to the first image by applying the OECF with 
the at least one parameter to the second image. 

♦ 

2. The method of claim 1, further comprising, prior to said generating a first histogram and 
generating a second histogram: 

removing a percentage of the overlapping pixels with the greatest difference in 
brightness. 

3. The method of claim 1, wherein: 

said generating a first histogram comprises recording in a first plurahty of pixel 
value bins a first plurality of numbers of pixels that have respective pixel values 
in the first region; and 

said generating a second histogram comprises recording in a second plurality of 
pixel value bins a second plurality of numbers of pixels that have the respective 
pixel values in the second region. 



4. The method of claim 3, wherein said detemiining corresponding pixel values in the first 
and the second histograms comprises generating a lookup table (LUT) storing a third plurality of 
nimibers of pixels and their corresponding pixel values. 

5. The method of claim 4, wherein said generating a lookup table coinprises: 

(1) initializing all entries in the LUT to 0; 

(2) initializing a first loop by setting / = 0;^ = 0; reml = hl[0]; and rem2 = h2[0]; 

, (3) updating the LUT by setting minj-em = mm{reml, rem2); reml = reml - 
min_rem\ and rem2 = rem2 - minj-em'^ and incrementing LUT[z][/] by minj-em; 

(4) if reml = 0, then incrementing / and setting reml = hl[i]; * 

(5) if rem2 = 0, then incrementing^ and setting rem2 = h2[/]; 

(6) if / < 256 md j < 256, then repeating steps (3) to (5); 

wherein hi [] is the number of pixels having a certain pixel value in the first 
histogram, h2[] is the number of pixels having a certain pixel value in the second 
histogram, and LUT[][] is the number of pairs of corresponding pixel values 
having a certain pixel value in the first histogram and a certain pixel value in the . 
second histogram. 

6. The method of claim 5, wherein said determining at least one parameter of an OECF 
comprises minimizing a color matching error, the color matching .error being defined as: 

255 255 

e = Y2LUTfiJfjJ((i + l)/256.0-S-'(TS((j + J)/256.0))), 

wherein e is the color matching error, r is a color matching parameter, and S() is the OECF. 

7. The method of claim 6, wherein the OECF is defined as: 
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wherein x is a pixel value normalized to (0,1), and a € (-1J ) is another color matching parameter. 

8. The method of claim 7, wherein said minimizing a color matching error comprises 
performing a golden section search of the color matching error. 

9. The method of claim 8, wherein said performing a golden section search comprises: 

(1) initializing the another color matching parameter and a minimum color 
matching error; 

(2) performing the golden section search with the another color matching 
parameter being fixed and the color matching parameter being varied to determine 
a smallest color matching error achieved; 

(3) recording values of the color matching parameter and the another color 
matching parameter that achieve the smallest color matching error if it is less than 
the minimimi color matching error; , 

(4) setting the minimum color matching error equal to the smallest color matching 
error; and 

(4) repeating steps (2) to (4) for a range of values of the another color matching 
parameter. 

10. The method of claim 6, wherein said applying the optoelectronic conversion function 
comprises: 

x,=S-^(W(t\xJS(xJ). 

wherein Xq is an original pixel value in the second image, Xc is a corrected pixel value in the 
second image, S^Q is the inverse of the OECF, and PTis a weight funciton defined as: 

W(r,x,)^T+(l-T)x,, 

1 1 . The method of claim 1 , wherein the OECF is defined as: 



wherein SQ is the OECF, jc is a pixel value normalized to (0,1), and a € (-1J j is a first color 
matching parameter. 

12. The method of claim 1 1 , wherein said determining at least one parameter of an OECF 
comprises minimizing a color matching error defined as: 

wherein e is the color matching error, xj and X2 are corresponding pixel values in the first and the 
second regions, Rj and i?2 are the first and the second regions, SQ is the OECF, S(X^ is the 
inverse OECF, and r is a second color matching parameter. 

13. The method of claim 12, wherein said minimizing a color matching error comprises 
performing a golden section search of the color matching error. 

14. The method of claim 13, wherein said performing a golden section search comprises: 

(1) initializing the first color matching parameter and a minimum color matching 
error; 

(2) performing the golden section search with the first color matching parameter 
being fixed and the second color matching parameter being varied to determine a 
smallest color matching error achieved; 

(3) recording values of the first and the second color matching parameters that 
achieve the smallest color matching error if it is less than the minimum color 
matching error; 

(4) setting the minimum color matching error equal to the smallest color matching 
error; and 
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(4) repeating steps (2) to (4) for a range of values of the first color matching 
parameter. 

15. The method of claim 12, wherein said applying the OECF comprises: 

x,=S-^(W(T.xJS(xJ), 

wherein Xo is an original, pixel value in the second image, Xc is a corrected pixel value in the 
second image, and Wis a. weight function defined as: 

W(t,xJ = t+(1-t)x,, 

16. A method for color matching a first image and a second image, wherein a first region of 
the first image and a second region of the second image overlap, the method comprising: 

removing a percentage of overlapping pixels with the greatest difference in 
brightness; 

generating a first histogram of the first region and a second histogram of the 
second region after said removing; 

histogram matching the first and the second histogram to determine corresponding 
pixel values from the first and the second histograms; 

minimizing a color matching error between the corresponding pixel values, 
wherein the color matching error is generated fi-om an optoelectronic conversion 
function (OECF); and 

color matching the second image to the first image by applying the OECF to the 
second image. 

17. The method of claim 16, wherein said histogram matching the first and the second 
histograms comprises generating a lookup table (LUT) as follows: 

(1) initializing all entries in the LUT to 0; 
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(2) initializing a first loop by setting i = 0;y = 0; reml = hl[0]; and rem2 = h2[0]; 

(3) updating the LUT by setting minj-em = mmfjemU rem2)\ reml = reml - 
min_rem; and rem2 = rem2 - minjrem; and incrementing LUT[/][/] by min_rem\ 

(4) if re/n7 = 0, then incrementing / and setting reml = hl[i]; 

(5) if rem2 = 0, then incrementing^ and setting rem2 = h2[/]; 

(6) if i < 256 and j < 256, then repeating steps (3) to (5); 

wherein hi [] is the number of pixels having a certain pixel value in the first 
histogram, h2[] is the number of pixels having a certain pixel value in the second 
histogram, and LUT[][] is the number of pairs of corresponding pixel values 
having a certain pixel value in the first histogram and a certain pixel value in the 
second histogram. 

18. The method of claim 17, wherein said determining parameters of the OECF comprises 
minimizing a color matching error, the color matching error being defined as: 

255 255 ... 

e = Y^^LUTfiJfJJ((i + 1)/ 256.0- S-'(rS((j + 1)/ 256.0))) , 

1=0 y=o 

wherein e is the color matching error, r is a color matching parameter and S() is the OECF 
defined as: 



2 

S(x) = x-^—arctan 
n 
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wherein SQ is the OECF, jc is a pixel value normalized to (0,1), and a e ) is another color 
matching parameter. 

19. The method of claim 1 8, wherein said minimizing a color matching error comprises 
perforaiing a golden section search of the color matching error comprising: 



- 14- 



(1) initializing the another color matching parameter and a minimum color 
matching error; 

(2) performing the golden section search with the another color matching 
parameter being fixed and the color matching parameter being varied to determine 
a smallest color matching error achieved; 

(3) recording values of the color matching parameter and the another color 
matching parameter that achieve the smallest color matching error if it is less than 
the minimum color matching error; 

(4) setting the minimum color matching error equal to the smallest color matching 
error; and 

(5) repeating steps (2) to (4) for a range of values of the another color matching 
parameter. 

20. The method of claim 19, wherein said applying the optoelectronic conversion function 
comprises: 

x,^S-'(W(T,xJS(xJ). 

wherein is an original pixel value in the second image, Xc is a corrected pixel value of the 
second image, and PTis a weight function defined as: 

W(t ,Xj) = T'\-(1-t)Xj. 
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